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Abstract 

We consider the most general primordial cosmological perturbation in a uni- 
verse filled with photons, baryons, neutrinos, and a hypothetical cold dark 
matter (CDM) component within the framework of linearized perturbation 
theory. We present a careful discussion of the different allowed modes, distin- 
guishing modes which are regular at early times, singular at early times, or 
pure gauge. As well as the familiar growing and decaying adiabatic modes and 
the baryonic and CDM isocurvature modes, we identify two neutrino isocur- 
vature modes. In the first, the ratio of neutrinos to photons varies spatially 
but the net density perturbation vanishes. In the second the photon-baryon 
plasma and the neutrino fluid have a spatially varying relative bulk velocity, 
balanced so that the net momentum density vanishes. Possible mechanisms 
which could generate the two neutrino isocurvature modes are discussed. If 
one allows the most general regular primordial perturbation, all quadratic 
correlators of observables such as the microwave background anisotropy and 
matter perturbations are completely determined by a 5 x 5, real, symmetric 
matrix-valued function of co-moving wavenumber. In a companion paper we 
examine prospects for detecting or constraining the amplitudes of the most 
general allowed regular perturbations using present and future CMB data. 

I. INTRODUCTION 



A key challenge of modern cosmology is understanding the nature of the primordial 
fluctuations that eventually led to the formation of large scale structure in our universe. 
One possibility is that the fluctuations were generated in a period of inflation prior to the 
radiation dominated era of the hot big bang. As inflation ended the fluctuations would 
then have been imprinted as initial conditions for the radiation era on scales far beyond the 
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Hubble radius. The second possibility is that the structure was generated through some 
causal mechanism operating within the standard big bang radiation and matter eras. In 
this paper we focus on the first option, that the fluctuations were imprinted early in the 
radiation era as linear fluctuations in the metric and in the matter and radiation content. 

For several good reasons the possibility that the primordial perturbations were adiabatic 
has been the focus of most interest to date. If the relative abundances of different particle 
species were determined directly from the Lagrangian describing local physics, one would 
expect those abundance ratios to be spatially constant because all regions of the universe 
would share an identical early history, independent of the long wavelength perturbations. 
The stress-energy present in the universe would then be characterized on large scales by a 
single, spatially uniform equation of state. Such fluctuations are termed adiabatic and are 
the simplest possibility for perturbing the matter content and the geometry of the universe. 
They are also naturally predicted by the simplest inflationary models fl]], although there also 
exist more complex models giving other types of perturbations [0. For a recent discussion 
see ref. 

Nevertheless, there is no a priori reason why the situation could not be more compli- 
cated with abundance ratios varying from place to place. Perturbations of this sort are 
known as isocurvature, or sometimes entropy, perturbations. Most studies of isocurvature 
perturbations have examined the possibility that the primordial perturbations were entirely 
isocurvature with a vanishing primordial adiabatic component and have sought to explore 
whether such pure isocurvature models could explain the observed structure in the universe 

m-m- 

In this paper we adopt a more phenomenological approach, which we believe is now war- 
ranted by the prospect of upcoming precision measurements of the CMB anisotropy on small 
scales. Ground based and balloon borne telescopes and the MAP and PLANCK satellites 
will provide very detailed measurements of the primordial fluctuations Most work on 

how to interpret this data has focused on parameter estimation starting from the assump- 
tion that the initial perturbations were generated by an inflationary model specified by a 
small number of undetermined free parameters 0. Adiabatic perturbations, characterized 
by an amplitude and spectral index, as well as tensor fluctuations also characterized by two 
parameters are usually assumed, and based on these assumptions a host of cosmological 
parameters, such as fltotai, ^b, ^a, h, N^, are to be inferred from the observed CMB multi- 
pole moments. While the prospects of such measurements are beguiling, they rely heavily 
on assumptions regarding the form of the primordial perturbations. We feel that those 
assumptions are worth checking against the data using an approach that assumes neither 
inflation nor any other favorite theoretical model. 

For learning about the fundamental physics responsible for structure formation, deter- 
mining whether the primordial fluctuations were in fact Gaussian and adiabatic is at least as 
important and interesting as measuring the values of cosmological parameters. For this pur- 
pose the relevant question is not whether primordial isocurvature perturbations offer a viable 
alternative to adiabatic perturbations, which has been the focus of prior work, but rather 
how and to what extent can observations constrain the presence of isocurvature modes. 
Rather than considering a particular isocurvature model, it seems appropriate to consider 
the most general primordial perturbation possible without adding new physics to the hot 
big bang era. In other words, we assume a universe filled with neutrinos, photons, leptons 
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and baryons, and a cold dark matter (CDM) component and then try to place constraints 
on the amplitudes of all possible perturbation modes. 

As mentioned above, we shall need to make a simplifying assumption to limit the param- 
eter space of possible perturbations to reasonable size. We shall assume that the fluctuations 
were indeed primordial — that is, that the perturbation modes were excited by physics op- 
erating at a very early epoch preceding the hot radiation dominated era and that no new 
dynamics influenced the perturbations at late times. Here 'late' means well before decou- 
pling, so that whatever decaying modes were excited earlier had a chance to decay before 
influencing the observed structure in our cosmic microwave sky. This assumption excludes 
cosmic defect models (i.e., strings, global textures, etc.) |jl2|] of structure formation where a 
detailed understanding of the dynamics of the order parameter field is required to determine 
the unequal time correlations at late times p!^JI^ . 

How should the most general perturbation be characterised? This question has his- 
torically generated some confusion in the literature. The standard terminology refers to 
'growing,' 'decaying,' and 'gauge' modes, the latter referring to modes affected by general 
coordinate transformations preserving the chosen gauge condition. In this paper as in most 
work on cosmological perturbations we shall use synchronous gauge. The variables in other 
common gauges (e.g., Newtonian gauge) are merely linear combinations of the synchronous 
gauge variables and their time derivatives. In synchronous gauge, the two gauge modes for 
scalar perturbations are easily identified. The first corresponds to time-independent spatial 
reparametrizations of the constant time hypersurfaces. In this mode the metric perturbation 
is constant and the matter is unperturbed. The second corresponds to deformations of these 
hypersurfaces through a spatially dependent shift in the time direction. For the latter gauge 
mode, which diverges at early times, both the matter variables and the metric are perturbed. 
Indeed this mode may be regarded as a shift in the time of the big bang singularity. The 
remaining modes are then defined modulo the two gauge modes. 

We characterise the remaining modes as either 'regular' or 'singular,' according to their 
behavior as the time since the big bang tends to zero. This terminology is preferable to 
the standard one because it includes constant modes such as the neutrino velocity mode 
we shall discuss. In this paper we shall treat only the 'regular' modes (i.e., those regular 
up to the gauge modes). There are several reasons for this. Singular modes are necessarily 
decaying as one proceeds forward in time away from the big bang, so if they are present at 
some very early time with linear amplitude so that a perturbative treatment is valid they 
quickly become irrelevant. A second reason is that even if the perturbation amplitude is 
small the perfect fluid approximation breaks down at early times for the singular modes. 
Higher moments in the Boltzmann hierarchy become progressively more important as one 
goes back in time, and specifying the initial conditions involves specifying an infinite number 
of constants. The perfect fluid approximation seems essential to a simple specification of 
primordial initial perturbations. Of course, decaying modes might be produced by some 
late time physics (such as cosmic defects) operating well after the big bang, but our point 
here is that it would appear very difficult to characterise them without introducing explicit 
source terms. In contrast, the 'regular' modes are completely characterised by specifying 
the leading terms in a power series expansion in conformal time of the low moments of the 
phase space density — that is, the density and velocity of the fluids. 

Having characterised the regular perturbations by the leading terms in the power series 
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for the metric and fluid densities and velocities, any quadratic observable (e.g., the matter 
power spectrum or the cosmic microwave anisotropy power spectrum) is then completely 
determined by a primordial power spectral matrix, which rather than a single function of k 
as it would be for growing mode adiabatic perturbations is a 5 x 5, real, symmetric matrix 
function of k. The off-diagonal elements establish correlations between the modes. As long 
as only quadratic observables are considered, no assumption of Gaussianity is required. 

The counting arises as follows. Each cosmological fluid is described by two flrst order 
equations, so each fluid introduces two new perturbation modes. However, in synchronous 
gauge, as mentioned, there is one gauge mode affecting the fluid perturbations. For a single 
fluid there is therefore just one regular, growing mode perturbation, and no physical (i.e., 
non-gauge) decaying mode. A convenient way to deal with the gauge mode is to identify 
it with the velocity of the cold dark matter. By a coordinate choice this may be chosen 
initially to be zero. If it is assumed that the cold dark matter couples to the other fluids 
only via gravity, there is no scattering term to consider, and with coordinates chosen so that 
the velocity is initially zero, it remains so for all times. If we now introduce photons and 
baryons, four new modes arise. The flrst is an adiabatic decaying mode. There are also the 
baryon isocurvature and cold dark matter isocurvature modes, where the initial conditions 
contain equal and opposite perturbations of the photon density and the baryon or cold dark 
matter densities, respectively. The fourth new mode is that in which the photon and baryon 
fluids have a relative velocity that diverges at early times and the fluid approximation breaks 
down. So far we have three regular modes. Let us now introduce neutrinos into the picture. 
We shall imagine we are setting up the perturbations after neutrino decoupling (T ~MeV, 
t ~ seconds) . For our purposes there is no distinction between the different neutrino species 
nor between neutrinos and antineutrinos, since we are only interested in how perturbations in 
the neutrino fluid affect cosmological observations of the density and microwave background 
today. Two new perturbation modes are introduced, the flrst a neutrino isocurvature density 
perturbation and the second a neutrino isocurvature velocity perturbation. In the latter, we 
can arrange the neutrino and photon-baryon fluids to have equal and opposite momentum 
density. In the approximation that we ignore the collision term coupling neutrinos to the 
photon-baryon fluid, which is valid after neutrino decoupling, and in the small time limit 
there is no divergence in this mode. 

The neutrino isocurvature velocity perturbation may be considered as primordial as long 
as one takes primordial to mean 'generated after one second' but well before photon-baryon 
decoupling. Of course similar remarks apply to the cold dark matter, baryon isocurvature 
or neutrino isocurvature perturbations. Namely if we go back far enough in cosmic history, 
where the various conservation laws for baryon or lepton number break down, or where the 
cold dark matter was initially generated or reached thermal equilibrium, then a description 
of the form we are using would also break down. Since one second is still quite early, and 
certainly well before photon-baryon decoupling, we regard it admissible to consider the pos- 
sibility of 'primordial' neutrino velocity perturbations. Possible mechanisms for generating 
them shall be discussed in a separate paper P2| . 

A primeval baryon isocurvature model (FBI) was introduced by Feebles [0 in which a 
universe with just baryons, radiation, and neutrinos is assumed and primordial perturbations 
in the baryon-to-photon ratio are assumed. Since at early times the baryons contribute 
negligibly to the total density, such perturbations lack an adiabatic, or curvature, component 
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at early times. Compared to an ri = 1 CDM model PBI gives: (1) lower small-scale 
relative peculiar velocities, (2) greater large-scale flow velocities, (3) earlier reionization, 
and (4) earlier galaxy and star formation 0. The consequences of PBI and comparison 
with observations have been studied by a number of authors . 



Bond and Efstathiou [|T3,§ have considered a CDM isocurvature model in which the 
CDM-to-photon ratio varies spatially. A possible mechanism for generating such pertur- 
bations arises in models with axion dark matter in which scale invariant fluctuations of 
the axion field A{x) are converted into density fluctuations after the axion field acquires 
a mass in the QCD phase transition |]T3|]. Under the assumption that quantum fluctua- 
tions during inflation impart a scale free spectrum of fluctuations to the axion field, so that 
5A{k) ~ fc"^/^, and that SA{k) <^ A, a scale-free spectrum of Gaussian fluctuations in the 
axion-to-photon ratio is generated. The resulting spectrum of density fluctuations today has 
the same power law on large scales as for adiabatic fluctuations with a scale free spectrum 
Pp{k) ~ k^, but compared to adiabatic, scale-free perturbations the turnover to Pp{k) ~ k~^ 
power law behavior on small scales occurs on a larger scale for the isocurvature case. For the 
amplitude of density perturbations normalized on small scales (for example, using ag), this 
gives about thirty times more power in the matter perturbations on large scales and entails 
an excessive CMB anisotropy on large angles. This work has been extended to consider 
mixtures of CDM isocurvature and adiabatic fluctuations, low density universes, and more 



recent CMB data p|,|T7| 



The neutrino isocurvature modes discussed here allow for a spatially varying relative 
density of photons and neutrinos and a relative velocity between the photon and neutrino 
components as well. For the neutrino isocurvature density mode, the total density pertur- 
bation vanishes but the relative density of neutrinos and photons varies spatially. On super- 
horizon scales, the neutrinos and photons evolve similarly but upon entering the horizon the 
neutrinos free stream, developing nonvanishing higher moments of the neutrino phase space 
density Fyi{k), while because of Thomson scattering the photons continue to behave much 
as a perfect fluid. These distinct behaviors subsequently generate perturbations in the total 
density. For the neutrino isocurvature velocity mode, the rest frames of the neutrino and the 
photon fluids do not coincide. The relative velocities are such that initially the perturbation 
in the total momentum density vanishes. If this last condition were not satisfied, the metric 
perturbation generated by this mode would diverge at early times, rendering the mode a 
singular mode, which according to the discussion above should be excluded. But the lack of 
divergence owing to this cancellation makes this an admissible regular perturbation mode. 
The neutrino isocurvature mode solutions are implicit in the work of Rebhan and Schwarz 
P^l and of Challinor and Lasenby but their implications were not explored. 

A possible obstacle to exciting neutrino isocurvature modes arises at early times, from 
processes in which neutrinos are generated or scattered. If the neutrino chemical potentials 
vanish, so that there are in each generation precisely as many neutrinos as antineutrinos, 
a spatially varying relative density can only be established at a temperature sufficiently 
low that the processes turning photons into neutrino-antineutrino pairs and vice versa had 
already been frozen out — that is, below a few MeV. Nonvanishing chemical potentials for 
the various neutrino species can protect variations in the ratio {p^/ p^) from erasure by pro- 
cesses involving uu annihilation, and observational constraints that would rule out neutrino 
chemical potentials of the required magnitude are lacking. L-violating processes mediated 
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through sphalerons, unsuppressed at temperatures above the electroweak phase transition, 
can readjust the neutrino chemical potentials /ij,^, fi^^, and /ly^ and, moreover, can con- 
vert lepton asymmetries in the neutrino sector into baryon asymmetry, but if the neutrino 
chemical potentials satisfy fi^^ + + = 0, then no tendency favoring sphalerons over 
anti-sphalerons is introduced and the net effect of these electroweak processes vanishes. Re- 
call that the neutrino overdensity is proportional to fii^J + ii^^^ + /i^y^^. Similarly, for the 
neutrino isocurvature velocity mode scattering of neutrinos with other components dampens 
this mode at very early times, and for this mode to be relevant there must exist a process 
capable of exciting it after this dampening effect has frozen out. 

There exist many candidate mechanisms that might generate these neutrino isocurvature 
modes. The neutrino density isocurvature mode could be generated during inflation if the 
theory included a light scalar field carrying lepton number, with mass much smaller than 
the Hubble constant during inflation, During inflation this field would be excited but 
would contribute negligibly to the density of the universe. After inflation, when the Hubble 
constant fell below the mass of the scalar field it would oscillate and decay, producing a lepton 
asymmetry. We would expect the neutrino chemical potential to be proportional to the scalar 
field value, so the most natural possibility would be Gaussian, scale invariant perturbations 
in the neutrino chemical potential, with the density perturbation in the neutrinos being 
proportional to the square of the chemical potential. 

The neutrino velocity mode could have been excited by the decay of relics such as cosmic 
strings, walls, or superstrings, surviving until after the neutrinos had decoupled and then 
decaying into neutrinos. The neutrino fluid produced in such processes would be perturbed 
both in its density and velocity, and these perturbations would be isocurvature in character 
||22|| . Another possibility for exciting the modes arises from magnetic fields frozen into the 
plasma whose stress gradients impart a velocity to the photon-lepton-baryon plasma but 
not to the neutrino component. 



II. IDENTIFYING THE MODES 



We now turn to identifying the possible perturbation modes using synchronous gauge 
with the line element 



ds = a 



dr^ + {5ij + hij) dx^dx^ 



For spatial dependence of a given wavenumber k, we define 



hij{k, r) 



fc.^,Mk,r)+(A:.A:,--5.,)6,(k,r) 



Our discussion generally follows the notation of ref. ||TT 
The linearized Einstein equations are 



(1) 



(2) 



k^T] - -nh = (AnGa^) 5T%, 

h + 2nh - 2k'^r] = -{SnGa^) 6T\, 
h + 6ii + 27i{h + 6?)) - 2k'^r] = -{24:7iGa'^){p + p)a 



(3) 
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where we define {p + p)9 = ik^ T^j and (p + p)cr = —(kikj — ^5ij)Ti^ . We define Ti = a/a. 
With only a single fluid, 6 is simply the divergence (V ■ v). 

Assuming N fluids, labeled ( J = 1, . . . , A^), one may rewrite the above as 

k'v-lnh = ~n'Y.^j6j, 

h + 2nh - 2k'^r] = -m^ J2 ^J, 

J 

h + 6ii + 2nih + 6?)) - 2k'^ri = -UTfVt^a^ (4) 

where wj = pj/ pj and = dpj/dpj. In the last equation only the neutrino contribution to 
the anisotropic stress is included. At early times only the neutrino contribution is relevant, 
but later when the photons and baryons begin to decouple the photon quadrupole moment 
must also be included. 

For the photons, the equations of motion at early times before the baryonic component 
of the fluid was significant are 



4 2 ■ 

3 ^ 3 



5^ + -0^ + -/i = 0, 



- = 0. (5) 

Similarly, for the neutrinos, after neutrino decoupling at temperatures of ~ 1 MeV, we have 

4 2 . 

+ -Ou +3^ = 0' 

9u - \kH, + Pa, = 0, 
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2d, + h + Qr]\-^kF,^. (6) 



where = Fy2/2 is the quadrupole moment of the neutrino phase space density and F^i is 
the /-th multipole, in the notation used below. 

The photon and neutrino evolution equations above differ only in the presence of an 
anisotropic stress term a^. Photons, because of their frequent scattering by charged leptons 
and baryons, at early times are unable to develop a quadrupole moment in their velocity 
distribution. Neutrinos, on the other hand, develop significant anisotropic stresses upon 
crossing the horizon. The addition of an extra degree of freedom that would result from the 
equation for is avoided by setting 0",^ = at r = 0.Q 

For a CDM component, the equations of motion are 



If we were to consider the physical decaying mode with r] oc t , we would find this condi- 
tion cannot be satisfied. As discussed above, decaying modes are inevitably associated with the 
breakdown of the fluid approximation at early times. 
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Sc + dc + ^h = 0, 

e, + ne, = o. (7) 

At later times, when Qfy becomes comparable to Q^, the photon equations of motion are 
modified as follows: 

4 2 • 

5^ + + 3^ = 0' 

- —6^ + an^aT^O^ - Oh) = 0. (8) 

Here a is the scale factor, Ue is the electron density, and ax is the Thomson scattering cross 
section. The baryon equations of motion are 

Sb + Ob + ^h = 0, 
0b + nOb + ^^aneCTTiOb - e^) = 0. (9) 

At early times, the characteristic time for the synchronization of the photon and baryon 
velocities tf,^ ~ l/(ne(JT) is small compared to the expansion time t^xp ~ and to the 
oscillation period tosc ~ (ar/k). (We use units where the speed of light c = 1.) Any 
deviation of (^^^ — 6b) from zero rapidly decays away. This may be seen by subtracting the 
second equation of ^ from that of ^ and regarding TC9b + as a forcing term. In the 

limit (Tt — > cxD, one obtains 9b = 6^ — in other words, a tight coupling between the baryons 
and the photons. Therefore, at early times we may set 9b — ^7 — 9^b- 

In the tight coupling approximation the evolution equation for 9^b is obtained by adding 
the second equation of | multiplied by to the second equation of |^ multiplied by 
so that the scattering terms cancel giving 

+ 9^b = -^bn9^b + \n^k'^5^ (10) 
and the baryon and photon density contrasts evolve according to 

'5b = —9^b — -h, 

6, = -^9,b - h. (11) 

In the absence of baryon isocurvature perturbations, 6^ = making one of these equations 
redundant, but both equations are required for the most general type of perturbation. 

Although an excellent approximation early on, the tight coupling assumption breaks 
down at later times as the photons and baryons decouple, and for a more accurate descrip- 
tion the two equations in ^ for the time derivatives of the monopole and dipole moments of 
the velocity of the photon phase space distribution must be replaced by the following infi- 
nite hierarchy of equations for the time derivatives of higher order moments of the photon 
distribution function as well: 
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3 ' 



0, = k'{^-^)+an,aT{e,-e,), 



3, „ 4 ; 8 . 9 

k 



7^ 



2i + l 



- + l)-F7(^+i) - aUeaxF^e (12) 



where £ > 3 for the last equation. Initially, as r — > 0, F^i = for £ > 2. If this condition 
were relaxed an infinite number of decaying modes would appear. The fluid approximation 
for the photons, eqn. H, is obtained using only the first two of the above equations combined 
with the approximation F^2 = 0. 

To describe the neutrinos during and after horizon crossing requires a Boltzman hierarchy 
for 6^, 6^, F^i identical to the one above except that the Thomson scattering terms are 
omitted. This assumes that neutrino masses are irrelevant. 

Before solving the Einstein equations ^, we first identify the two gauge modes result- 
ing from the residual gauge freedom remaining within synchronous gauge. A general in- 
finitesimal gauge transformation considered to linear order is a coordinate transformation 
x'^(x) where = x'^ + e^(x). From the transformation rule for tensors 

9'^..i^') = ^^3^^'^, (13) 

it follows that 

r £ ^3^} . £ (0) , (0) 5 

= + e\M glu + (14) 

where g^^J is the unperturbed, zeroth-order metric. Applying ^ to the metric ds"^ = 
(^'^{'t) —(1 — ^oo)'^T^ + S/iQj drdx^ + {6ij + hij)dx^dx^ , where 

t = t' + T{^,t), x = x' + S(x,r) (15) 

where S and T are regarded as infinitesimal, one obtains 

5^00 = -—T{^, r) - 2T(x, r), 
a 

Shoi = -rj(x, r) + Si{x, r), 

6hij = Sij + Sj^i H 6ijT. (16) 

For the density contrast for a perfect fluid component labeled by a, one obtains S{Sa) = 
—3(1 + Wa)'HT. After a gauge transformation the velocity is shifted by (5v = — S'(x, r) for 
all components. 

The synchronous gauge condition 6hoo = implies that T(x, r) = A(x)/a(r), and 6hoi = 
implies that 
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S(x,r) =5(x) + VA(x) I 



dr 



Therefore 



dr 



air 



(17) 



(18) 



For a radiation-dominated universe [with a(r) = r] considering for the moment only a 
single wavenumber, we find that^ 



Shi 



6ij - 2kikj ln(r) A{k) + 2B kikj 



(19) 



or 



h = A 



- 2^ ln(r) 



2B, ri = A 



(20) 



To obtain a more accurate power series expansion close to the time of matter-radiation 
equality and to consider the baryon and CDM isocurvature modes, we assume a scale factor 
evolution for a universe filled with matter and radiation: a(r) = r + r^. At matter-radiation 
equality a^g = 1/4 and r^g = {\/2 — l)/2. Since H = (2r+ 1)/(t^ + t) has a pole at r = — 1, 
the resulting power series solutions are expected to diverge beyond the unit circle. 



For the matter-radiation universe, with a(r) = r + eqn. |20| is modified to become 



A 



r2 1 + r 



1 + r 



2B, 



V 



-A 



(l + 2r) 
r2(l + r)^ 



(21) 



We now present the power series solutions. and Ri, represent the fractional contribu- 
tion of photons and neutrinos to the total density at early times deep within the radiation 
dominated epoch. (We ignore possible effects due to nonvanishing neutrino masses). We 
also assume that we are considering the perturbations after the annihilation of electrons and 
positrons, so that the latter have dumped their energy into the photon background. For 
Ni, species of massless neutrinos we define R = |A^i/(^)^ and we have R^ = {1 + R)'^ and 
Ru = R{l + R)-\ 



Adiabatic Growing Mode. 



2^2 



h = -k'r 



7] = 1 



h + AR, 



2^2 



12(15 + 4/?,) 



Note that the gauge mode A disagrees with eqns. [94] and [95] of ref. [11 1, which are incorrect. 
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36 



1 

36 



23 + AR^ 



3(12 + 



15 + AR^ 



(22) 



Baryon Isocurvature Mode. 



h = AQbfiT - 6QbflT^, 
2 2 

5c = -2Qb,0T + 3^6,0^"^, 

^6 = 1 - 2f]6,oT + 3^l6,oT^ 
g 

g 



^76 — "o^&.O^^T^, 



There is no regular baryon velocity mode because of the tight coupling of the baryons to 
the photons. 



CDM Isocurvature Mode. 





4Oc,0T - 6Oc,0T^, 


V = 


2 2 




1 - 2Qc,oT + 31]c,or' 




—2Qc,oT + 3nc,oT^, 


= 


g 
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^c = 0, 

1 

_( 

3' 



^76 — ~77^c,0^^7"^, 



The CDM velocity mode may be identified with the gauge mode. This can be seen from the 
equation for the CDM velocity which is decoupled from the other equations. 6c dm behaves 
as at early times. If there were several species of CDM, however, a new physical, 
non-gauge relative velocity mode would arise, which would be divergent at early times. 



Neutrino Isocurvature Density Mode. 

h — -—k T 

R,, 



6(15 + 47?,, 



5e = ^^^k\\ 



207?^ 



r 1 Ri' 7 2 2 

Ofe = -—kr, 
8 R^ 

S ~ + —^^k^T^ 
R^ 6 7?^ 

5. = 1 - I^t', 
6 

a, = ^ ^ -A;V. (25) 

2(15 + 47?,) ^ ' 

Physically, one starts with a uniform energy density, with the sum of the neutrino and 
photon densities unperturbed. When a mode enters the horizon, the photons behave as a 
perfect fluid, while the neutrinos free stream, creating nommiformity in the energy density, 
pressure, and momentum density, thus sourcing metric perturbations. 



Neutrino Isocurvature Velocity Mode. 
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^~ 3(5 + 4i?,) ^ V (5 + 4i?,)(15 + 4i?,)y 

3r2bo-R!^(-R7 + 2) 2 

Oft = KT ■ KT , 

5^ = -:i^ifcr - ^bfiRujR^ + 2) ^_2 



3 -R^ -R2 

4 , ^b,oRu J 2 



= -i^kr ^ — /cr , 

^c = 0, 

Ru , , •i^bflRi/ , , -Ri/ /Sf^b o 9f2ft o^\ , 2 , R-i^ 1 3 2 



-f^'-y I^-y \^ J 



^^"^ 6(5 + 4i?.)^^' 

a = ^ kr + 16-^- 

3(5 + 4i?,) ^ (5 + 4i?,)(15 + 4i?,) ' 

^ fcVl (26) 



u3 



7(5 + 4i?,) 



Here the neutrinos and photons start with uniform total density and uniform density ratio 
but with relative velocities matched so that initially the total momentum density vanishes. If 
the relative momenta were not perfectly matched, the metric perturbation generated would 
diverge at early times, as in the adiabatic decaying mode. But because of the perfect match, 
a divergence at early times is avoided. 

It is also possible, at least in principle, to consider regular modes with higher moments 



of F^i with i > 3 excited initially, as was considered in ref. [25|; however, it is difficult to 



envision plausible mechanisms for exciting these higher moment modes. 

Newtonian Potentials. In this paper we have used synchronous gauge, but for complete- 
ness we give the form of the Newtonian potentials for the regular modes presented above. 
The Newtonian potentials and ip are related to the synchronous variables as follows: 

(P = ^[h + 6v + n(h + 6v)], 

^ = V-^[h + H- (27) 

We define the Newtonian potentials according to the convention ds^ = a2(r)[— (ir2(l + 20) + 
dx^ dx^ 6ij{l — 2ip)]. We now give the Newtonian potentials to leading order. For the growing 
adiabatic mode 

(15 + 4i?,)' 

^ = 7 7- (28) 

^ (15 + 4i?,) ^ ^ 
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For the neutrino isocurvature density mode 



For the neutrino isocurvature velocity mode 

-4:Ru 



_ —2Ry 
~ (15 + 4^ ' 

^ = TTF^^- (29) 



(15 + 4i?,; 



(15 + 4/2, 

The potentials for the baryon isocurvature mode are 

^ {AR, - 15)r]b,o 



^ = 7^^^.k-'r''- (30) 



-r, 



-r. 



^ 2(15 + 2/2,) 

. Z(i5.zi^,, (31) 
^ 6(15 + 2/?,) ' ^ ^ 

and for the CDM isocurvature mode are 

(4/2,-15)^, 
2(15 + 2/2,) 

It is curious that the Newtonian potential diverges at early times for the neutrino isocur- 
vature velocity mode while in synchronous gauge there is no singularity. It appears that 
synchronous gauge is a more physical gauge and that Newtonian gauge is inadequate for 
modes based on anisotropic stresses. The dimensionless Ricci curvature o?{t)t'^R is non- 
singular at early times. In any case, synchronous gauge is more physical in the sense that 
its evolution is local whereas in Newtonian gauge the evolution of the shape of the constant 
cosmic time hypersurfaces depends on how matter behaves infinitely far away (because the 
gauge choice is defined in terms of the scalar-vector-tensor decomposition, which is nonlo- 
cal). The divergence of the Newtonian potentials arises from the need to warp the surfaces 
of constant cosmic time so as to put the metric in a spatially isotropic form; however, there 
is nothing at all physical about this particular form. The neutrino isocurvature velocity 
mode was excluded in ref. because of the behavior of the Newtonian potentials at early 



times; however, when a physical phenomemon can be described in a nonsingular manner in 
some gauge, its singularity in another gauge should be regarded as a coordinate singularity. 



III. DISCUSSION 

We have identified five nondecaying modes corresponding to each wavenumber: an adi- 
abatic growing mode, a baryon isocurvature mode, a CDM isocurvature mode, a neutrino 
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10 100 1000 

FIG. 1. CMB Anisotropy for the Neutrino Isocurvature Density Mode. We plot 
lit + l)Q/27r for the neutrino isocurvature density mode (the dashed curves) for initial power 
spectra P^^ ~ A;" where a = —3.0, . . . , —2.4 increasing in increments of 0.1 from bottom to top at 
large t. The adiabatic growing mode (the solid curve) with a scale-invariant spectrum is included 
for comparison. All curves are normalized to COBE. For the lowest curve the variations in the 
photon-to-neutrino ratio obey a scale-invariant initial power spectrum. 



3x10"" - 
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J 

FIG. 2. CMB Anisotropy for the Neutrino Isocurvature Velocity Mode. We plot 
1(1 -t- l)Q/27r for the neutrino isocurvature velocity mode (the dashed curves) for initial power 
spectra ~ A;° with a = —3.0, . . . , —2.0 increasing in increments of 0.2 from bottom to top at 
large £. a = —2.0 corresponds to a white noise initial power spectrum in the divergence of the 
velocity field, possibly resulting from a large number of small explosions. 



15 



density isocurvature mode, and a neutrino velocity isocurvature mode. Under the assump- 
tion that the primordial perturbations are small enough so that the linear theory suffices, 
two-point derived observables are completely determined by the generalization of the power 
spectrum given by the correlation matrix 

{A^{k)A^{-k)) (33) 

where the indices (m, n = 1, . . . , 5) label the modes, independently of whether or not the 
primordial fluctuations were Gaussian. In Figs. 1 and 2 the shape of the CMB moments 
for the neutrino isocurvature modes arc indicated. In this computation the cosmological 
parameters Hq = 50km s~^, il^ = 0.05, Vic = 0.95 were assumed. In a companion paper, 
we examine prospects for constraining the amplitudes of these modes using upcoming MAP 
and PLANCK data. 
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